namespace algo {
    db E(Sta &sta) { // O(n^2)
        db E1 = 0, E2 = 0;
        for (Circ &x : sta) E1 += protDist(x);
        for (auto it1 = sta.begin(); it1 != sta.end(); ++it1)
            for (auto it2 = next(it1); it2 != sta.end(); ++it2)
                E2 += intrDist(*it1, *it2);
        return E1 + E2*E2;
    }
}
